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Abstract 

It has long been held that the low photosynthetic rates [A) of coffee leaves are largely associated with diffusive constraints 
to photosynthesis. However, the relative limitations of the stomata and mesophyll to the overall diffusional constraints to 
photosynthesis, as well as the coordination of leaf hydraulics with photosynthetic limitations, remain to be fully elucidated 
in coffee. Whether the low actual A under ambient CO2 concentrations is associated with the kinetic properties of Rubisco 
and high (photo)respiration rates also remains elusive. Here, we provide a holistic analysis to understand the causes 
associated with low A by measuring a variety of key anatomical/hydraulic and photosynthetic traits in sun- and shade- 
grown coffee plants. We demonstrate that leaf hydraulic architecture imposes a major constraint on the maximisation of the 
photosynthetic gas exchange of coffee leaves. Regardless of the light treatments, A was mainly limited by stomatal factors 
followed by similar limitations associated with the mesophyll and biochemical constraints. No evidence of an inefficient 
Rubisco was found; rather, we propose that coffee Rubisco is well tuned for operating at low chloroplastic CO2 
concentrations. Finally, we contend that large diffusive resistance should lead to large CO2 drawdown from the intercellular 
airspaces to the sites of carboxylation, thus favouring the occurrence of relatively high photorespiration rates, which 
ultimately leads to further limitations to A. 
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Introduction 

Because CO2 influx and water vapour efflux share a common 
pathway through the stomatal pores on leaf surfaces, a trade-olT 
between transpirational costs and CO2 assimilation is implicitly 
unavoidable. The couphng between stomatal conductance (gg) to 
CO2 and water vapour (and the need to maintain a proper leaf 
water balance) has often been evidenced by the strong positive 
scaling between and the leaf hydraulic conductance per unit 
area, Kj^ [1^3]. In turn, the significance of as a potentially 
limiting component of the vascular system has been further 
emphasised by the strong hydraulic-photosynthetic coordination 
observed across a large sample of diverse species [4] . Furthermore, 
is closely related to the anatomy of the leaf: Al has been shown 
to be positively related to both the theoretical axial conductivity of 
the midrib (determined from xylem conduit numbers and 
dimensions) and the venation density, [3,5]. A unified control 
of hydraulic and photosynthetic traits may also be further 
highhghted by comparing shade and sun leaves: the former have 



lower rates of net CO2 assimilation (A) and^^s, therefore leading to 
lower demand for water and correspondingly lower A'l and 
[4,6]. 

In addition to stomatal limitations, A is currently known to be 
constrained by mesophyll conductance {g,n), which is defined as the 
conductance for the transfer of CO2 from the intercellular 
airspaces (Q) to the sites of carboxylation in the chloroplastic 
stroma (Cy [7]. According to Flexas et al. [8], limitations to 
photosynthesis are of similar magnitude as stomatal constraints 
and generally greater than biochemical limitations. Increasing 
evidence has shown that ^n, is often intrinsically co-regulated with 
gs [8]. More recentiy, Ferrio et al. [9] showed a positive scaling 
between g^ and Al and proposed that water and CO2 share an 
important portion of their respective diffusion pathways through 
the mesophyll; thus, any downregulation of leaf hydraulics may 
reduce not only g^ but also ^ni, both of which contribute to 
reducing A. Quantification of g^ has become important in 
predicting leaf photosynthetic parameters using the Farquhar- 
von Caemmerer-Berry (FvCB) model of leaf photosynthesis [10] 
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because such a model underestimates the maximum Rubisco 
carboxylation rate (Kmax) by considering as being infinite 
[11,12]. Furthermore, changes in Q -Q because of variations in 
jgni among species result in different A values in plants with the 
same biochemical activity and gs [7] . 

At lower C„ A is limited by RuBP carboxylase activity, which in 
turn depends on the concentration, activation and kinetic 
properties of Rubisco [13]. For example, Rubisco with a higher 
specificity factor, 6',./„ (which determines the relative rates of 
carboxylation and oxygenation by Rubisco at given COj and O2 
concentrations), could be regarded as conferring an advantage in 
minimising photorespiration rates, Rp [14]. Interestingly, Rubisco 
seems to have evolved towards higher Sc/o under stressful 
conditions leading to low Q, such as drought, salinity and high 
temperature [14]. Additionally, species evolved under stressful 
conditions with sclerophyll leaves tend to display a higher 5, /,,, 
which may again be related to lower C[. [14,15]. However, trade- 
offs between S^/o (or particularly the affinity for CO2 (1/Ac)) and 
the carboxylase turnover rate (kc^) have been described, so that 
there is no Rubisco in nature with both high affinity for CO2 and 
fast activity [16]. At saturating Q., A becomes limited by the 
capacity for the regeneration of RuBP (often dominated by the 
electron transport capacity) but also by S^/„ [17]. In any case, 
when studying the ecophysiology of plants, Rubisco kinetic 
parameters should not be considered a constant but rather as an 
active part determining the biochemical potential of plants to 
assimilate COj under optimal and suboptimal conditions. 

Coffee is an evergreen perennial shrub with hypostomatous 
leaves that evolved in the African forest understoreys and is thus 
considered a shade-demanding species. However, in many 
situations, modem coffee cultivars grow well without shade and 
even out-yield shaded coffee [18,19]. At atmospheric CO2 
concentrations and saturating light, coffee displays low A, typically 
in the range of 4—1 1 |j,mol CO2 m ^ s ' [20,21], which is in the 
lower range recorded for trees [22]. However, the maximum A 
values obtained in common A/Q curves can surpass 20 ^mol CO2 
m~^ s~' [23], whereas the photosynthetic capacity, determined 
under true CO2 saturation (~50 mmol CO2 mol~' air), reaches 
values exceeding 30 |J,mol O2 m ^ s ' [21]. Taking the above 
information together, the low A in coffee might largely result from 
diffusive constraints to photosynthesis [24]. However, the relative 
contributions of the stomata and mesophyll to the overall 
diffusional limitations to photosynthesis, as well as the coordina- 
tion of leaf hydraulics with photosynthetic limitations, remain to 
be fully elucidated in coffee. Furthermore, the anticipated low 
and g^ would compromise the transference of CO2 through 
stomata and leaf mesophyll to Rubisco sites, therefore decreasing 
Cc and ultimately favouring RuBP oxygenation in relation to 
carboxylation [25]. Under these conditions, Rubisco kinetic 
properties with a high affmity for GO2 may be critical to secure 
optimal carbon balance and yield. 

Although the g(;n(Tal descriptors of the ecophysiolog)' of sun- 
shade plants are well-known, an integrative analysis of these 
descriptors has rarely been performed in conjunction (hydraulics, 
stomata, mesophyll, Rubisco, etc). Here, we provide a holistic 
analysis to understand the causes associated with low A in coffee by 
measuring a variet)' of key anatomical/hydraulic and photosyn- 
thetic traits, including D^, Kj^, the actual and maximum theoretical 
gs fewmax) and several parameters from A/Q and A/ Q curves. We 
have centred our attention on coffee given that, in addition to 
being a highly important commodity [18], it could also be 
considered as a model plant for other important crops which 
evolved as understorey trees, such as cacao, citrus and tea; these 
crops are traditionally considered to have low A (seldom above 



10 nmol CO2 m~^ s~', even in the field under favourable 
conditions) likely as consequence of large diffusi\'e, rather than 
biochemical, limitations to photosynthesis [26]. Our main goals 
were (i) to examine the role of leaf hydraulics as a prime limiting 
factor for photosynthetic gas exchange, (ii) to calculate g^ to 
properly parami^terise th(; responses of A to C,., and (iii) to 
disentangle the relative contributions of stomatal, mesophyll and 
biochemical limitations to photosynthesis, and how all of these 
facts may be affected by the light supply. A fiirther goal was to 
analyse how the kinetic properties of coffee's Rubisco [27] may 
affect the actual A and evaluate if Rubisco in coffee is well tuned 
for operating at low C,. by comparing it with the highly efficient 
Rubisco of Limonium gibertii; this shade-intolerant species is 
particularly adapted to highly stressing environments, whose 
Rubisco has evolved under extremely low Cc ranges [14]. We 
demonstrate that leaf hydraulic architecture imposes a major 
constraint on the maximisation of the photosynthetic gas exchange 
of coffee leaves and that A was mainly limited by stomatal factors 
followed by similar limitations associated with mesophyll and 
biochemical constraints, regardless of the light supply. We also 
found that, despite the relatively high S^/o of coffee Rubisco, the 
large diffusive resistances favour the occurrence of relatively high 
Rp, which ultimately leads to further limitations to A. 

Materials and Methods 

Plant Material, Growth Conditions and Experimental 

Design 

The experiment was conducted in Vi?osa (20°45'S, 42°54'W, 
altitude 650 m) in southeastern Brazil. Uniform seedlings of Cojfea 
arabica L. cv 'Catuai Vermelho lAC 44' obtained from seeds were 
grown in 12 L pots containing a mixture of soil, sand and 
composted manure (4:1:1, v/v/v). Plants were grown either under 
fuU sunlight conditions (100% light) or under low light in a shade 
environment (10% full sunlight). The shade enclosure was 
constructed using neutral-density black nylon netting, and the 
plants were kept in these conditions for 12 months before 
measurements. Throughout the experiment, the plants were 
grown under naturally fluctuating conditions of temperature and 
relative air humidity and were fertilised and irrigated as necessary. 
The pots were randomised periodically to minimise any variation 
within each light environment. For all samplings and measure- 
ments, the youngest fuUy expanded leaves, corresponding to the 
third or fourth pair from the apex of the plagiotropic branches, 
were used. The experiment was arranged in a completely 
randomised design, with six plants in individual pots per treatment 
as replicates. The experimental plot included one plant per 
container. 

Morpho-anatomical Features and Related Hydraulic and 

Diffusive Traits 

The specific leaf area (SLA) was computed using the dry mass of 
20 leaf discs (1.13 cm diameter each). For anatomical measure- 
ments, leaves were collected and fixed in FAA70 for 48 h, followed 
by storage in 70% (v/v) aqueous ethanol. Samples of the mean 
region of each leaf blade were embedded in methacrylate 
(Historesin-Leica Microsystems Nussloch. rh-idcllx-rg, Germany) 
according to the manufacturer. Transverse sections (7 )j.m 
thicknc-ss), obtained using a rotary microtome (model RM2155, 
Leica Microsystems Inc., Deerfield, USA), were stained with 
toluidine blue at pH 4.0 and mounted in synthetic resin 
(Permount). Anatomical data were quantified using an image 
analysis program (Image Pro-Plus, version 4.5, Media Cybernetics, 
Silver Spring, USA). Video images were acquired using a video 
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camera attached to an Olympus microscope (AXyOTRT, Olym- 
pus Optical, Tokyo, Japan). The following anatomical data were 
then assessed: (i) total leaf thickness; (ii) palisade and spongy 
mesophyll thickness; (iii) upper and lower epidermis thickness; (iv) 
guard cell length (Z); (v) vertical distance from the vein to stomatal 
epidermis, D^.^:, (vi) stomatal density (SD, according to DaMatta 
et al. [28]); (vii) stomatal index that was estimated as 
(S/(E-|-S)) X 100, where S and E are the number of stomata 
and epidermal cells per unit leaf area, respectively [29]; and (viii) 
stomatal pore index based on guard cell length (SPIgei) that was 
estimated as SPIgci = SD x Z,^. For determining venation traits, 
samples of leaf lamina were cut from two leaves per plant and 
cleared as described by Scolfoni et al. [30]. Regions of approx- 
imately 6 mm^ were imaged at SOX, and was calculated as the 
sum of the vein lengths divided by the total image area using the 
image analysis program described above. Leaf size was estimated, 
using the maximum leaf widths and lengths, with the equations 
described by Antunes et al. [3 1] . Stomatal size and epidermal cell 
size, as well as the level of coordination between anatomical traits 
and leaf size (quantified as the deviation from the expected 
proportional relationship to each other), were all determined 
following Carins Murphy et al. [6]. 

Midrib xylem conduits were measured to determine the 
theoretical midrib axial hydraulic conductance {K^, where the 
conduits were treated as ellipses to calculate Ki as 

where a and b are the long and short internal vessel diameters, and 
ri is the viscosity of water at 25°C [32], further normalising by leaf 
length and area. 

The maximum stomatal conductance to water vapour (g,vmax) 
was calculated according to Franks et al. [33] as 

max = SD a/[y(J+n/2y/{a/n)] 

where SD is the stomatal density, rfw is the diffusivity of water 
vapour in air, a is the maximum area of the open stomatal pore, v 
is the molar volume of air, and / is the stomatal pore depth for fully 
open stomata. Values for standard constants rfw a^^d v were those 
for 25°C (24.9x10"'^ m^ s"' and 24.4x10"^ m^ mor\ 
respectively), a was calculated as TiipjTf' , where p is the stomatal 
pore length, which was approximated as LI 2 according to Franks 
and Farquhar [34]. / for fully open stomata was taken as Z/4, 
assuming guard cells inflate to a circular cross section [33]. 

Leaf Hydraulic Conductance 

The leaf hydrauhc conductance (A'l) was estimated according to 
Brodribb and Holbrook [35] by following the kinetics of water 
potential ("Pi) relaxation in rehydrated leave as: 

J^L = Cln(<i'o/>Pf)/f 

where C is leaf capacitance, estimated using pressure-volume 
curves [36], *Po is *Pi before rehydration, and *Pf is T] after 
rehydration for t seconds. Leaf Tj was measured using a 
Scholander-type pressure chamber (model 1000, PMS Instru- 
ments, Albany, NY, USA). 

Gas Exchange and Fluorescence Measurements 

Leaf gas exchange and chlorophyU a fluorescence were 
measured simultaneously using an open-flow infrared gas- 



exchange analyser system equipped with a leaf chamber fluorom- 
eter (LI-6400XT, Li-Cor, Lincoln, NE, USA). Environmental 
conditions in the leaf chamber consisted of a leaf-to-air vapour 
pressure deficit between 1.2 and 2.0 kPa and a leaf temperature of 

25°C. 

In light-adapted leaves, the actual quantum yield of PSII (<I>psii) 
was determined by measuring steady-state fluorescence (i^s) and 
maximum fluorescence during a light-saturating pulse of c. 
8,000 |imol m~^ s ' (Zin'), following the procedures of Genty 
et al. [37]: 

®PSn = (^m'-J's)/J'm' 

The electron transport rate (7f) was then calculated as 

7F = QiPPPFDa)psii 

where PPFD is the photosyntheticaUy active photon flux density, a 
is the leaf absorptance and P is the PSII optical cross section. The 
product a P was herein determined from the relationship between 
<l>PSii and <l>f;o2j obtained by varying the CO2 concentration 
under non-photorespiratory conditions in an atmosphere contain- 
ing less than 1 % O2, as described by Valentini et al. [38] . 

The light respiration rate (/?l) was determined according to the 
'Laisk-method' [39], using thej axis intersection of^/C; (internal 
CO2 concentration) curves performed at three different PPFD 
intensities (750, 250 and 75 |J,mol m~^ s~'). The rate of 
mitochondrial respiration at darkness {Rq) was measured early in 
the morning in dark-adapted leaves. 

The photorespiratory rate of Rubisco (i?p) was obtained 
according to Valentini et al. [38] using the following equation: 

iJp = l/12[/F-4(^ + -RL)] 

Six A/ C; curves were obtained from different plants per 
treatment. In light-adapted leaves, Al Q curves were initiated at 
an ambient CO2 concentration (CJ of 400 |a,mol moP' under a 
saturating PPFD of 1000 jtmol m"^ s"'. Once steady state was 
reached, C, was decreased stepwise to 50 |xmol moP' air. Upon 
completion of the measurements at low C„ was returned to 
400 nmol mor' air to restore the original A. Next, was 
increased stepwise to 2,000 jtmol mor' air. Al Q curves consisted 
of 13 diflerent values. 

Cc was calculated after Harley et al. [40] as 

Co = (r* (7f + 8(^ + Rl)))/{Jf - A{A + 

where Y* was determined from the in vitro Rubisco specificity 
factor (i'c/o) as 

r* = 0/5c/o 

A was taken from gas-exchange measurements, and the J-g 
values were obtained from chlorophyll a fluorescence yield. After 
estimating Q, was calculated following Harley et al. [40] : 
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g„ = ^/(Ci-Co) 

From the A/ Cj and Al curves, the maximum carboxylation 
capacity (F^max) and maximum capacity for electron transport rate 
(7miix) were calculated on a Q and Q. basis using the kinetic 
parameters for coffee described in Martins et al. [27]. The FvCB 
model was fitted to the data by applying iterative cur\ c- fitting 
(minimum least square difference) using the Microsoft Excel Solver 
tool (Microsoft Corporation, Redmond, WA, USA). Additionally, 
&i) ^cmax and Jmax were estimated using the Ethier and Livingston 
[1 1] method, which is based on fitting A/ Q curves with a non- 
rectangular hyperbola version of the FvCB model, relying on the 
hypothesis that gm reduces the curvature of the Rubisco-limited 
portion of an ^/C; response curve. Again, the kinetic parameters 
of Rubisco measured on colfee were used. Corrections for the 
leakage of COg and water vapour into and out of the leaf chamber 
of the Li-6400-40 have been applied to all gas-exchange data, as 
described by Rodeghiero et al. [41]. 

The chloroplastic CO2 concentration of transition (Q ,rans); 
where Q denotes the transition between the Rubisco- and RuBP 
regeneration-limited states, was estimated as described by Gu 
et al. [42]: 

^c_trans ~(«^max-^m ^^cmaxr )/(4F^i3iax "Anax) 

where is the effective Michaelis-Menten constant for CO2 that 
considers the competitive inhibition by Oj which was taken from 
Martins et al. [27]. 

The overall photosynthetic limitations were partitioned into 
their functional components [stomatal (4), mesophyU (4J and 
biochemical (4>) limitations] using the values of^s, gm, P^maxj T , 
and Cc following the approach proposed by Grassi and Magnani 
[43]: 

4 = feot/gsX02 * 5^/5Cc)/feot + dA/dC^ ) 



Im = (gtot/gm * dA/dCo)/(gtot + dA/dCo) 



/b =gtot/(ftot + SA/ dCc) 

where co2 is thfi stomatal conductance to CO2 
{gs_C02=gs/i-6), is the mesophyU diffusion conductance 
according to Harley et al. [40] and ^„„ is the total [:onductance 
to CO2 from ambient air to chloroplasts 
(ftot = l/[(l/^s-co2) + (l/gm)]). dA/8Cc was calculated as: 



dA/dCo = [Vc max(r* +K,n)]/iCc + K,^f 



Statistical Analyses 

Data are expressed as the means ± standard error. Student's t- 
tests were used to compare the parameters between treatments. 
Additionally, one sample <-tests were performed to compare the 

means for shade plants with the expected values if they were 
proportional to the sun plants. AU statistical analyses were carried 
out using Microsoft Excel. 



Results 

Sun- and shade-grown individuals differed reasonably in 
venation architecture and mesophyU structure (Table 1). The 
sun leaves, compared with shade leaves, displayed a higher leaf 
thickness (15%) primarUy resulting from a higher thickness of 
palisade (43%) and spongy (14%) mesophyU, which led to a higher 
palisade-to-spongy parenchyma ratio (25%) and lower SLA (39%). 
However, the upper and lower epidermis thickness and the 
maximum distance from the vein to the epidermis (i)v-e) did not 
change in response to the light treatments. Sun-grown plants had 
31% more midrib xylem conduits but with slightly lower mean 
conduit diameters (7%) than shade-grown plants. The sun plants 
displayed and Kj^ values that were 160% and 58% higher, 
respectively, than in the shade individuals. Adjustments to the light 
treatments associated with the morphological characteristics of 
stomata were also observed, as denoted by higher SD (62%), 
stomatal index (27%)) and SPIg,.| (62%) in sun leaves compared 
with the shade leaves; however, both the guard cell length and 
stomatal size did not change in response to the light treatments. 
Vein density was proportional to 1/V leaf area but slighdy higher 
(17%) than would be expected if it was directiy proportional to SD 
(Figure 1). In line with the differences in the stomatal index, SD 
and epidermal ceU size deviated significantiy from the expected 
proportional relationship to leaf area: SD was 25%) higher and 
epidermal ceU size was 3 1 % lower in shade plants than would be 
expected if these variables were directiy proportional (Figure 1). 
The maximum g^ (g„max)> determined by the stomatal size and 
density, was larger (61%) in sun leaves than in the shade leaves. 

Under saturating PPFD (1,000 |tmol photons m~^ s~'). both A 
and gs values were higher (c. j5"/o) in sun leaves than in the shade 
leaves (Table 2). Tlu' values for several other photosynthetic/ 
respiration parameters were also higher in sun leaves than in the 
shade leaves (Table 2): 60% for g^ (average of two independent 
mediods), 42% for Fc„ax, 45% for J^^ (for V.^^.^^ and J^^, tiie 
values express averages obtained on Q and Q bases), 140% for 
R-D, 200% for ^L, 41% for R^, and 51% for Jf. In contrast, C\, C„ 
Ccjrans ^ud the ratios A/g^, g^/gs and Jmax/ Vcma^ on both Q and 
Q bases did not differ significantiy in response to the light 
treatments. Fcmax and Jmax calculated on a Q basis were, on 
average, 101% and 37% higher than P^max andj^ax! respectively, 
calculated on a Q basis, whereas the Jmax/^cmax ratio was on 
average 32% lower when calculated on a Q basis than on a Q 
basis. 

The overall photosynthetic limitations were essentially simUar 
when comparing sun- and shade-grown plants (Table 2). We 
found that A was mainly constrained by stomatal limitations (c. 
40%)), whereas mesophyU and biochemical limitations accounted 
simUarly for the remaining limitations {c. 30% each). 

The analysis of Rubisco kinetic properties are summarised in 
Table 3, showing that coffee presents a Rubisco with a relatively 
high affinity for CO2 (low K^j and fast activity (ICcm'^, resulting in 
relatively high values for Sc/o- From the combination of in vivo 
^crnax and in vitro data {kcJ^, the concentration of active Rubisco 
sites was estimated to be 16 and 20 jtmol m~^ s~' for shade and 
sun plants, respectively. 

Discussion 

This study provided a hohstic examination of the key steps that 
could limit the photosynthetic capacity to fix CO2 and demon- 
strated that leaf hydraulic architecture imposes a major constraint 
on the maximisation of the photosynthetic gas exchange of coffee 
leaves by limiting It is tempting to suggest that these constraints 
might to some extent explain the photosynthetic behaviour of 
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Figure 1. The relationships between vein and stomatal density, vein density and 1/v leaf area, stomatal density and 1/leaf area, and 
epidermal cell area and leaf area. Data (± standard error) are shown for coffee plants grown under shade (blacl< circles) or full sun (white circles). 
Asterisks denote differences (* P<0.05; **, P<0.01) between the observed values for the shade plants and the expected values (broken line) if these 
were proportional to the averages for the sun plants. 
doi:10.1371/journal.pone.0095571.g001 



other important (sub)tropical crops such as cacao, citrus and tea, 
which have a photosynthetic performance and a slow-growth 
behaviour similar to that of coffee [26] . Our results suggest that 
improvements of the photosynthetic performance of these crops by 
means of selecting hydraulic traits (e.g., Z)„ and Kj) that might 
support higher values (thereby decreasing stomatal limitations of 
photosynthesis) could be a useful strategy to facilitate the selection 
of promising genotypes with enhanced crop growth and produc- 
tion. 

We have shown that adjustments in leaf hydraulics through 
increases in and in sun-grown individuals were coordinated 
with a suite of traits related to water flux and gas exchange per leaf 
area, such as mesophyU structure (e.g., higher mesophyU thickness 
and palisade-to-spongy parenchyma ratio) and stomatal attributes 
(increased SD and stomatal pore index). Importantly, we showed 
that increases in Kt occurred at the expense of an increased 
number of midrib conduits with lower lumen in sun plants, 
suggesting that improvements in the hydraulic safety would not 
compromise the hydraulic efficiency [44]. In addition, the larger 
SD coupled with larger stomatal index implies a proportionally 
greater increase in the number of guard cells than in normal 
epidermal cells [45]; hence, light availability should directly and 
positively influence stomatal fate in coffee regardless of passive 
changes linked by light-induced differences in leaf blade expan- 
sion. However, differential leaf expansion was responsible for the 
adjustment of to SD, which remained nearly proportional to 



each other. Thus, coffee plants can balance water supply with 
transpirational demand through a coordination of increased 
initiation of stomata cells with differential expansion of epidermal 
cells. Such coordination reflects an optimisation of the trade-off 
between transpirational costs and CO2 assimilation, resulting in 
the higher intrinsic water use efficiency observed in coffee [c. 
85 ixmol CO2 moF' H2O on average), which is markedly high 
compared with other tropical woody species [46,47] . 

Our A'l values, determined using the method of relaxation 
kinetics of ^i, were remarkably lower than those found by 
Brodribb et al. [1] for other tropical trees using the same method 
(between 17 and 36 mmol m ^ s ' MPa '). Considerably low 
values of Al had already been reported for C. arabica [c. 4.2 mmol 
m ^ s ' MPa ') by Gasco et al. [48] using the high-pressure 
method. Consistent with the close relationship between ATl and 
[4,5,49], our maximum Z)„ values were also within the low range 
recorded for angiosperms [50], suggesting that the hydraulic 
architecture of coffee leaves imposes strong resistance for water 
flow which, in turn, should limit the CO2 diffusion into the leaf 
[4] . In addition, our ATl and D^, values (Table 1) fit very well in the 
general relationship of both A"l and with the maximum A 
proposed by Brodribb et al. [5,51], that is, our estimated Al and 
Dv values exactiy predict our maximum measured A values. 
Therefore we contend that the leaf hydraulic architecture 
ultimately should act as a prime factor limiting photosynthetic 
gas exchange in coffee leaves. 
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Table 1. Anatomical and hydraulic traits of coffee plants grown under shade or full sunlight conditions. 



Parameters 


Treatments 






Shade 


Full sunlight 


Specific leaf area {m^ i<g"^) 


22.9±3.5 


14.0±2.7* 


Total leaf thickness (fim) 


333.9±9.1 


384.5 ±17.2* 


Palisade thickness (^im) 


52.6±1.2 


75.24±4.4* 


Spongy thickness (jim) 


220.9 ±6.1 


253.1 ±11.5* 


Upper epidermis thickness {\im) 


38.8±0.2 


37.4±1.3 


Lower epidermis thickness (nm) 


29.6±0.9 


26.4 ±1.2 


PP/SP 


0.24±0.01 


0.30±0.01* 


SPIgd 


0.104±0.006 


0.168±0.003* 


Guard cell length (|im) 


28.0 ±0.7 


28.9±0.4 


Stomata! density (mm^) 


129.1 ±7.2 


208.8±4.2* 


Stomatal index 


20.3 ±1.1 


25.7±1.0* 


gwmax (mol m"^ s"') 


1.16±0.05 


1.82 ±0.02* 


Venation density (mm mm"^} 


4.0±0.1 


5.5±0.2* 


Midrib vessel diameter (^m) 


22.8 ±0.2 


2 1.3 ±0.4* 


Number of midrib conduits 


117±7 


153±11* 


Kt (mmol MPa"' m"^ s"') 


24.5 ±2.5 


63.8±1.6* 


Kl (mmol MPa"' m"^' s"') 


6.9±0.1 


10.9 ±1.4* 


Dv-e (tim) 


195.8±8.6 


215.4±8.5 



n = 6±SE. Asterisks indicate statistically significant differences (P<0.05) between shade and full sun treatments. 

Abbreviations: PP:SP, palisade-to-spongy parenchyma ratio; SPIgd, stomatal pore index based on guard cell length; gwmax, maximal theoretical stomatal conductance to 
water vapour; Kt, midrib conductance; Kl, leaf hydraulic conductance; Dv-e, vertical distance from the vein to the stomatal epidermis. 
doi:l 0.1 371 /journal.pone.0095571 .tOOl 



Regardless of light regimens, our actual values were relatively 
similar to the value of 108 mmol H2O m s averaged over a 
range of studies using cofiee plants grown under optimal 
conditions ([10], and references therein). These values are 
significantiy lower than those of the modelled g^^,^, resulting in 
a gy^i^/gs ratio above 13, as can be calculated for sun-grown 
plants. By comparison, that ratio was c. 2.5 in non-water-stressed 
tomato [52] and 4.5 in Eucalyptus globulus [33]. The high hydraulic 
resistance of the coffee leaves is most likely the cause of the 
difference between the theoretical g„j„.^ and the recorded g^. 
Nevertheless, this large difference raises the following questions: 
why would the plant invest in a large ^vmiax if the maximum 
realisable A is relatively low and constrained by the leaf hydraulic 
architecture? Furthermore, from an ecological point of view, what 
would be the advantage of having a large g^msoi if, in the humid, 
shaded understoreys where coffee evolved, A should be more 
constrained by light limitation than by CO2 supply? Despite not 
having immediate responses to these queries, our results suggest a 
lack of coordination between the maximum capacity for stomatal 
aperture and carbon fixation, as also noted in saplings of Bornean 
rainforest tree species grown in the understorey [53]. A large gy^^^ 
may not be problematic in terms of water loss in the humid 
understorey, where transpiration rates are expected to be much 
more dependent on boundary layer resistance, and thus the 
importance of the stomata in optimising photosynthetic gas 
exchange should be reduced. In any case, considering that bigger 
stomata tend to close slower than smaller stomata [33,54], the 
relatively large stomata of coffee leaves (combined with low Ki) 
might result in excessive leaf desiccation if large stomatal apertures 
are realised. In this sense, the low actual g^ might be a conservative 
strategy to maintain leaf hydration and minimise the risk of xylem 



embolism. This observation is in line with recent results obtained 
for Toona ciliata, where transpirational homeostasis to changes in 
vapour pressure deficit was achieved through dynamic stomatal 
control rather than modification of the relationship between veins 
and stomata [55]. Taken together, these findings lead to the 
interesting question of why long-term adjustments to the 
parameters that define g^i-nax have not been frxed and regulation 
of gs comes predominantly from short-term adjustments to 
environmental conditions but at the cost of inherendy low g^ in 
some species, such as coffee. 

Our maximum g^-^ value was c. two-fold higher than those 
previously reported for C. arabica seedlings [56]. In any case, our 

values were similar to those obtained for other evergreen woody 
species (e.g. [57-59]). Greater g^ values for sun leaves, as found 
here, have been systematically reported and have been often 
explained by anatomical and morphological differences between 
shade and sun leaves [60]. Despite the changes in g^ and A 
between shade- and sun-grown coffee plants, C; and Q remained 
fairly similar. Thus, regardless of the light treatments, when A 
changed, and g^ scaled accordingly, and, hence, C,. remained 
constant [61]. This proportional scaling lends support to explain 
why stomatal, mesophyU and biochemical limitations to photo- 
synthesis were similar between sun and shade leaves. These 
findings are in agreement with other studies, which show that the 
approximate scaling of g^ and g^-^ with A makes the relative 
limitations to photosynthesis rather conservative between sun and 
shade leaves, as also noted in Fagus sylvatica [62] . 

Irrespective of the light environment, the mean drawdown from 
C; to C,- [c. 79 |J,mol mol ) was lower than that from to Q (c. 
153 limol mol~'), which is consistent with the fact that the 
diffusive limitations to COj in mesophyU were lower than those in 
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Table 2. Mean values for the photosynthetic 
conditions. 


and respiration parameters of coffee plants grown 


under shade or full sunlight 




Parameters 


Treatments 






Shade 


Full Sun 


A (\imo\ CO2 m"^ s"') 


7.9 ±0.4 


12.0±0.8* 


Qs {mmol H2O s~^) 


94±9 


146±17* 


gm_Harley (mmol CO2 m"^ S"') 


76±4 


116±10* 


gm_Ethier (mmol CO2 m"^ s"') 


65±12 


109±5* 


Cj ([imol CO2 mol"^ air} 


247+11 


246 ±7 


Cc_Hariey (jimol CO2 mol"' air) 


143+10 


142±9 


Jf (fimol rm"^ s"') 


70±5 


106±5* 


Vcmax„a (nmol CO2 m"^ s"') 


26.8 ±1.6 


42.1 ±2.2* 


l/cmax_Cc (jimol CO2 m"^ S"') 


58.1 ±3.5 


78.5 ±4.0* 


Vmax.ci (nmol e" CO2 m"^ s"') 


71.1 ±4.4 


109.5±3.1* 


-'max.cc (fimol e" m"^ s"') 


1C2.7±6.3 


142.6±7.4* 


-^max/i/cmax_a 


2.6 ±0.04 


2.7±0.1 


^ max/^cmax_Cc 


1.8±0.1 


i.8±a.i 


Q.trans (^moi CO2 mol"' air) 


228±17 


249 ±23 


Rd (nmol CO2 m"^ s"') 


0.5 ±0.06 


1.2 ±0.05* 


Rl (|imol CO2 m"^ s"' 


0.1 ±0.05 


0.3 ±0.08* 


/A/g, (|imol CO2 mol"' H2O) 


86±6 


84±4 


gm„Harley/gs (mol COj mol"' CO2 ) 


1.4±0.1 


1.3±0.1 


Rp (jimoi CO2 m"^ s"') 


3.2 ±0.4 


4.5±0.1* 


Stomatal limitation 


0.41 ±0.04 


0.38±0.02 


Mesophyll limitation 


o.3a±o.a2 


0.30±0.01 


Biochemical limitation 


0.29±0.03 


0.32±0.03 



n = 6±SE. The overall photosynthetic limitations associated with stomatal, mesophyll and biochemical factors are also shown. Data for A, g^, gm_Hariey' Q_Hariey' -^f- 3nd 
Rp were obtained under PPFD of 1000 jimol m"^ s"\ Ca of 400 (imol moP' and leaf temperature of 25 'C. Asterisks indicate statistically significant differences (P<0.05) 
between shade and full sun treatments. 

Abbreviations: A, net photosynthesis; g^, stomatal conductance to water vapour; gm, mesophyll conductance to CO2; C|, sub-stomatal CO2 concentration; Q, 
chloroplastic CO2 concentration; Jp- electron transport rate estimated by chlorophyll fluorescence; L^cmax- maximum carboxylation capacity; Jmax- maximum capacity for 
electron transport rate; Q trans- the Q that denotes the transition between the Rublsco- and RuBP regeneration-limited state; Rp, dark respiration; Rl, light respiration; Rp_ 
photorespiration rate. 
doi:l 0.1 371 /journal.pone.0095571 .t002 



stomata, thus ultimately reflecting a high ratio. Indeed, by 

analysing the dataset reported by Flexas et al. [63], we noted that 
species operating at low Q tended clearly to display an increased 
&n/& ratio (Figure SI) which contributes to an improved 
performance in photosynthetic water use [63]. We believe that 
such feature is essential for keeping a positive carbon balance given 
that stomatal limitations may even be exacerbated in coffee trees 
grown under field conditions, particularly because peaks in the 
early morning and decreases progressively throughout the day, 
reaching values typically ranging from 10 to 50 mmol H2O 



m ^ s ' from midday onwards [23,24,64—67] as a consequence of 
rising vapour pressure deficit. 

Under saturating PPFD, A at ambient CO2 was limited by 
Rubisco activity regardless of the light treatment given that the 
estimated Q was lower than Q trans (Table 2). However, the 
realised A and J^p at maximum growth irradiance of shade plants {c. 
200 |J.mol PPFD m ^ s ') are c. 60% of their saturating values 
(Figure S2), indicating that, even though these plants operate 
during most of their development under light limitation, the 
balance between RuBP regeneration and Rubisco activity (the 
Jmax/ K:majc ratio) was csscntially similar in the shade and sun 



Table 3. Rubisco kinetic constants measured for coffee (taken from Martins et al. [27]). 



Species 5^,^ T* (Hbar) (^lM) K„ (^llVl) 



Coffee 



n = 4±SE. 

Abbreviations: S^/o, Rubisco specificity factor; F*, CO2 compensation point In the absence of mitochondrial respiration; and K^, the Michaells-Menten kinetics for CO2 
and O2, respectively; k^^^, Rubisco catalytic turnover rate for the carboxylase reaction. 
doi:l 0.1 371/journal.pone.0095571 .ta03 
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relatively low Q., which could, to some extent, reduce the revenue 
stream in a scenario of increasing CO2 atmospheric levels. 

The rates of respiration and photorespiration are other key 
processes influencing the plant carbon balance. In sun plants, R^, 
and Rp corresponded to 10% and 38% of A, respectively. 
However, given that i?L represented only a fraction (25%) of 
which is in line with the inhibition of mitochondrial respiration in 
the presence of light [71,72], the realised constraint of respiration 
on A is expected to be relatively low. In turn, Rp is expected to 
significandy affect A at midday, when the stomatal closure in coffee 
leaves exacerbates the drawdown from Q to Q [24], thus 
favouring the oxygenase activity of Rubisco. 

Conclusion 

Regardless of the light treatments, A was mainly limited by 
stomatal factors followed by similar limitations associated with the 
mesophyll and biochemical constraints. Our data suggest that an 
increased gm/gs ratio coupled with a Rubisco well- tuned for 
operating at low Q might be adaptations to a lower Q resulting 
from the low g^, these adaptations might have helped the 
establishment of coffee plants when they were moved from the 
understoreys to the more stressful conditions of open fields, where 
high vapour pressure deficit and elevated temperatures may 
constrain gas exchange. Our results also suggest that the coffee 
plant does not properly optimise its resource allocation: on the one 
hand, there seems to be an over-investment in capacity for 
carboxylation and electron transport despite limited light avail- 
ability under shade; on the other hand, an excess of electron 
transport capacity is to be expected under full sun. Nevertheless, 
this excess might be useful in attempts to increase A in coffee 
through breeding aimed at improving hydraulic traits {D^ and Kj) 



35 




Figure 2. Comparison of tlie CO2 assimilation rate as a function of C„ modelled with the kinetic parameters of Rubisco in coffee and 
several C3 species. The RuBP saturated rates of CO2 assimilation at 25 "C were calculated using the following equation [13]: 
A = (Cc — T')-{kcaf-[Kuhisco])/{Cc + Kc{l + 0/Ko)} — R'L. Q, the chloroplastic CO2 partial pressure; 0, the intercellular O2 partial pressure; T*, 
the CO2 compensation point in the absence of day respiration; [Rubisco], the catalytic site content of Rubisco. The Rubisco kinetic properties for 
coffee (Table 3) and the kinetic properties for the other species were retrieved from Savir et al. [69] and are summarized in Table SI. The Rl and 
[Rubisco] were set as 1.0 |tmoi CO2 m"^ s"' and 25 nmol sites m"^ respectively. 
doi:1 0.1 371/journal.pone.0095571 .g002 



plants (Table 2). These data indicate that, in contrast to the 
optimal distribution principle, which suggests that Rubisco and the 
regeneration of RuBP should co-limit photosynthesis such that no 
excess capacities remain [68], the coffee plant does not properly 
optimise its resource allocation. On the one hand, under light 
saturation, an excess of electron transport capacity is to be 
expected, given that the actual Q operates far from Cc_trans (cf. 
Table 2). On the other hand, under light limitation, a great 
investment in capacity for carbon assimilation and electron 
transport is retained despite the low reahsed A and Jf by the 
shade leaves. 

Coffee Rubisco was characterised as displaying (i) a relatively 
high 6'c/o, similar to that of woody evergreen species from xeric 
habitats [14]; (ii) a higher affinity for CO2 (low A^.), which ranks 
coffee Rubisco with the third kiwest value oi R'^ among C3 and C4 
plants recorded to date [69]; and (iii) a relatively high k^..^t , 
superior to the average of species from warm climates [70] . We 
next modelled the responses of ^ as a function of Q by comparing 
these Rubisco properties with those of several C3 species (Table 
SI) including Limonium gihertii, the species with the highest S^./„ and 
lowest R]. reported to date, thus particularly adapted to low Q. 
[14]. Importantly, we found that, all else being equal, the A values 
that would be achieved using the Rubisco kinetic properties of 
coffee or Limonium would be quite similar and superior to those 
from all other species analysed (Figure 2) suggesting that coffee 
presents an "efficient" Rubisco well-tuned for operating at low C,.. 
Additionally, the higher kcm'' and lower Kc mean that fewer 
Rubisco molecules are required to realise a given A. Nevertheless, 
given that lower leads to a reduction in Q transi ^ of coffee 
leaves is expected to begin to be limited by RuBP regeneration at 
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and expectedly supporting a high g^. Finally, we contend that the 
large diffusive resistance should lead to large drawdown from C,^ to 
Cc, thus favouring the occurrence of relatively high despite the 
relatively high S^,/„ of coffee Rubisco, which ultimately leads to 
further limitations to A. 

Supporting Information 

Figure SI The relationship between sub-stomatal GO2 
concentration and mesophyll-to-stomatal conductance 
ratio in the multi-species dataset from Flexas et al. [63]. 

(TIF) 

Figure S2 The response of net photosynthetic rate and 
electron transport rate to the photosynthetic photon flux 
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